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Abstract 

There are some distinguished ensembles of non-Hermitian random matrices for 
which the joint PDF can be written down explicitly, is unchanged by rotations, 
and furthermore which have the property that the eigenvalues form a Pfafnan point 
process. For these ensembles, in which the elements of the matrices are either real, 
or real quaternion, the kernel of the Pfaffian is completely determined by certain 
skew orthogonal polynomials, which permit an expression in terms of averages over 
the characteristic polynomial, and the characteristic polynomial multiplied by the 
trace. We use Schur polynomial theory, knowledge of the value of a Schur polynomial 
averaged against real, and real quaternion Gaussian matrices, and the Selberg integral 
to evaluate these averages. 

1 Introduction 

Let X r be an N x N random matrix with elements independent and distributed as standard 
real Gaussians, and let X q be a 2N x 2N random matrix with successive two-by-two blocks 
of the form 

^jk ^jk 

with each Zjk and Wjk independently chosen standard complex Gaussians. Equivalently, 
the joint probability density function (PDF) of X r is proportional to e - TrXI ( xr ) T / 2 an d that 
of X q is proportional to e - TrXq ( xq ) t / 2 These are said to define the real and quaternion real 
Ginibre ensembles respectively [16] . Notice that the distribution of the former is unchanged 
by either of the mappings X r i— > X r O and X x i— > OX v , for O G O(N), while the latter 
is unchanged by either of the mappings X q i— >■ X q S and X q i— > SX q for S G Sp(2AT). 
Real and quaternion real matrix ensembles possessing these properties will be denoted 
ROtv(50 an d QS 2 7v(5') respectively, where g is the underlying joint PDF of the matrices. 
Thus R07v(e _TrXxT//2 ) corresponds to the real Ginibre ensemble, and QS 2 7v( e ~ TrXXt ^ 2 ) 
corresponds to the quaternion real Ginibre ensemble. 



1 



For both these ensembles, the corresponding eigenvalue distribution has the property 
that the corresponding fc-point correlation functions are given by a 2k x 2k Pfaffian, in which 
the elements depend on a certain kernel function independent of k. In the quaternion real 
case this was established in [20] and [TUl Exercises 15.9 Q.2], and was earlier conjectured by 
Mehta [25]. The understanding of this property for the real Ginibre ensemble came later, 
and is exhibited in the works [131 I2S1 12] . More recently ensembles of real and quaternion 
real matrices with joint PDFs proportional to 

(detgtg)^ 

det(I + GtG)°2 1 ' ' 

and 

(detG t G) /3l det(I-G t G) /32 , \\G*G\\ < 1 (1.2) 

have been shown to also have the property that their eigenvalues form a Pfaffian point 
process; see p2], [2H [7J in relation to the former, and [21], [I] for the latter. We remark 
that the joint PDF (11.11) . for certain a 1; a 2 gives the distribution of the matrix G = A^ X B, 
with A and B independent Ginibre matrices [22| [T2l 121"] . while the joint PDF (II. 2p . for 
certain Pi, p 2 gives the distribution of sub-blocks G of unitary matrices with real, complex 
and real quaternion elements [291 El EEH EH H] . It is also the case that these point processes 
can be naturally projected onto the sphere in the case of (II. ip [T8l fT2] and the anti-sphere 
in the case of (11.21) |llj . [23] . providing us with the names spherical and anti-spherical 
ensembles respectively. 

Another construction of random matrices with joint PDFs (II. ip and (ll.2p allows for 
G to be rectangular. In the case of (II. ip . let A = XX\ where X is an M x mi, mi > 
M standard Gaussian matrix with real (/3 = 1), complex (P = 2) or real quaternion 
(/? = 4) entries, and let Y similarly denote such Gaussian matrices, but now with mi 
replaced by N, N < M. Then we know from [17] and [10], Exercises 3.6 q.3] that the 
rectangular random matrix G := A~ X I 2 Y has joint PDF given by (11.11) with a,i — and 
«2 = (/5/2)(mi + N). In relation to (II. 2p the rectangular matrix is literally a rectangular 
sub-block of the unitary matrices with real, complex and real quaternion elements [291 Ej ; 
specifically, for a block G of size N x M from a K x K unitary matrix (when the elements 
are real quaternions, all these sizes should be doubled) the joint PDF (11.21) results with 
fa = and p 2 = (P/2)(K -N-M+l- 2/(3). Now form the square matrix U(GG^) 1 ' 2 , 
where U is a random unitary matrix with the elements real, complex or real quaternion 
as for G. This random matrix has the same joint distribution as G, but multiplied by the 
factor (det GG^) a with a = (M - N)/2 in the real case, and a = M — N in the complex 
and quaternion real cases [5] , and thus gives rise to (II. ip and (II. 2p with non-zero values of 
ai and f3i [6], [7J (see also Section 3.2 below; this class of square random matrices formed 
from rectangular random matrices is referred to as an induced ensemble ). 

An important point is that the generalized partition function for the joint eigenvalue 
PDF (in the real case this requires summing over the allowed number of real eigenvalues — 
0,2, ... ,N for iV even) can be expressed as a Pfaffian in which the matrix elements are 
a skew inner product (•, -) s of a polynomial basis of the monomials up to and including 
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x 2N ~ l ] in the real case see [MJ Eqns. 8-10], [I2J Prop. 3.1] and [2H Eqns. (3), (6)], while 
for the quaternion real case see [2D], [23], [7]. The skew orthogonal inner product for 
all the real cases of the above joint PDFs is the same except for the details of the one- 
body functions therein, and similarly for the quaternion real cases. Moreover, the kernel 
function for the correlation functions is completely determined once the polynomials are 
chosen as {Qj(x)}j = o,i,2,..., each polynomial Qj(x) monic and of degree j, and having the 
skew orthogonality property 

(Q2j,Q2k+l)s = r j3j,k (j > k) (Q2j,Q2k)s = (Q2j+1, Q2k+l) s = 0. (1.3) 

The common structure means that there is a common structure for the skew orthogonal 
polynomials. In the real case this was identified by Akemann et al. [U Eqns. (4.6)-(4.7)] 
and reads 

Qln( z ) = (det(^I 2n - G)) G6R o 2n 
Qln+ii*) = z Qln(z) + (det(*I 2n - G)TrG) GeRQ2n (1.4) 

(in fact in general Q2n+i( z ) is not unique — any multiple of Q2n( z ) can be added without 
affecting (II. 3p . Formally the same expressions hold in the quaternion real case, 

Ql n {z) = <det(*I 2n - G)) GeQS2n 
QLn(*0 = zQl£z) + (det(zl 2n - G)TrG) GeQS2 „, (1.5) 

a fact known earlier than for the real case, and due to Kanzieper |20j . 

It is our aim in this paper to show how to compute the matrix averages in (jl.4|) and 
(11. 5p for the real and quaternion Ginbre ensemble, its induced counterpart, as well as the 
real and quaternion real versions of the spherical and anti-spherical ensembles specified by 
the joint PDFs flU]) and (Ol) . 

2 Schur polynomials, matrix integrals and general- 
ized Selberg integrals 

Let G be a square matrix. The theory of group characters tells us that there a special 
symmetric polynomial in the eigenvalues of G — the Schur polynomial — indexed by a 
partition k := (k\ , . . . , kn) , and denoted s K (X). One relatively simple property of the 
Schur polynomials is the expansion 

TV 

det(zI N -G) = J2 z N - l (-l) l Sl i(G), (2.1) 

1=0 

where l' denotes the partition with 1 repeated / times and all other parts zero. This follows 
immediately once one is familiar with the fact that Sii(x) = ei(x), where x := (x\, . . . , x^) 
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and ei denotes the l-th elementary symmetric polynomial. The expansion ( 12. ip will be 
used to substitute for det(zl 2n — G) in ( II .4p and (jl.5p . 

This expansion in the second of the formulas in ( 11.40 and (jl.5p leads to the product 
of a Schur polynomial and Tr (X) =: e\(X). Such a product can be expanded in terms of 
Schur polynomials by appealing to the Pieri formula (see e.g. [TU1 Prop. 12.8.4]) 

e m {x)s K (x) = ^ s *( x )- ( 2 - 2 ) 

A 

A/k a vertical m strip 

Here X/k consists of those boxes of the diagram of A which are not in k and this is said to 
be a vertical m-strip if X/k consists of m boxes, all of which are in distinct rows. Thus 

Sl i(G)Ti (G) = s (2|ll -i)(G) + s ll+ i(G). (2.3) 



More sophisticated properties of the Schur polynomials are the evaluations of the matrix 
averages which result from these substitutions. In fact for X real and unchanged by 
X i-> XO and X ^ OX, O e O(N) we have [H Eq. (3.6)] (see [27] for the case A = I) 

M AX)) X AB^ {P " )(XXT))X ' ^ 
0, otherwise, 

where the partition Ik is the partition obtained by doubling each part of k, and Pi a \x) 
is the one parameter generalization of the Schur polynomial known as the symmetric Jack 
polynomial (the case a = 2 as appears in ( 12.41) can equivalently be defined as the zonal 
polynomials associated with the symmetric space gl(N,M)/0(N))\ see e.g. (101 Ch. 12]. 
Similarly for X quaternion real and unchanged by X y XS and X y SX S G Sp(2iV) 
we have O Eq. (3.8)] 

{ Pj 1/2) (AA T ) (1/2) t 2 
(pM\l*))* {K 1 /i " (2-5) 

0, otherwise, 

where k? is the partition obtained by repeating each part of k twice. 

The final main ingredient in our evaluations of the skew orthogonal polynomials are 
averages of Jack polynomials with respect to Selberg type densities. To present these 
formulas, let ME^jvQj^x)) denote a matrix ensemble with eigenvalue PDF proportional to 



N 

\Xk Xj I 



n»w n 

1=1 l<j<k<N 
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The Selberg weight corresponds to the case g(x) = x a (l — x) b , < x < 1. Furthermore, 
introduce the generalized Pochhammer symbol 



u 



T(n - (j -l)/a + Kj) 



=1 r(n- (j -l)/a) • W 



An integration formula conjectured by Macdonald and proved by Kaneko and Kadell (see 
PH Eq. (12.153)] tells us that 



fflsi a ,^)=Pi a) (0 [Al + (iV ~ 1)/a + ^TT - (2.7) 

2/Q ' Jv( ( } ; [ Al + A 2 + 2(iV-l)/a + 2]l a) 

Also required is the limiting case (A 2 — > oo and scaling of the integration variables) 

(PW(X!, . . • , x^jMB.^^e-.j = P^Cl^IAi + (N- I) /a + 1]W (2.8) 
and the variant [2H] (see also [3], [15] ) 

[ Xl + (N-l)/a+lp 

(-1)N[-A 2 + Ax + 2 + 2(N - I) /ap ' 

(2.9) 



(PW(X 15 . . . ,^)) ME2/Qjv( ^ 1(1+:c) -, 2) = PW(1^ 



3 Skew orthogonal polynomials for the Ginibre en- 
sembles 

3.1 The classical Ginibre ensembles 

Using the explicit form of the skew inner products, the skew orthogonal polynomials for 
the real, and quaternion real, Ginibre ensembles have been computed in earlier literature. 

Proposition 1. We have fWj . flfy Prop. 5.10.4] 

Ql n (z) = z 2n , Q\ n+l {z) = z 2n+1 - 2n^" 1 (3.1) 
and Jgg/, JM Exercises 15.9 Q.2] 

n 21 

= QUz) = 2 n n\ ^ (3-2) 

1=0 

Here we want to rederive these results by a direct evaluation of fll.4p and ( II. 5p . Consider 
first the real Ginibre ensemble. Substituting (12.1 p in the first formula of fll.4p gives 

In 

Qln(z) = Y, z2n ~ l (- i y( S ^ G »GeR0 2n{ e-^CCT /2 y (3.3) 
1=0 
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According to (12. 4p the only non-zero case of the average in (13.31) is / = 0, so giving 
Q v 2n (z) = z 2n as stated in f lXTj) . 

Substituting (12. ip in the second formula of (II .4p . then making use of (12. 3p tells us that 

(det{zl N - G)TrG) GeR02n{e _ TlGGT/2) 

2n 

= z 2n ~\-l) 1 : ((s (2)li -i ) (G f )) GeR02n{e ^ GG T /2) + (sii+i(G)) GeR02n(e _ TrGG T /2) y (3.4) 

1=0 

According to (12 .4p the only non-zero case of the averages in (13.41) is / = 1 in the first, and 
so we have 

(det(*Ijv - G)TrG) G£R02n{e ^ GG T /2) = (p* {GG T )) G( . R02n{e _^ GGT/2) . (3.5) 

Mi I 1 J 

Changing variables to H = GG T in the RHS of (13. 5p . which introduces the Jacobian 
factor (det if) -1 / 2 (see e.g. [TQl Prop. 3.2.7]), then changing variables to the eigenvalues and 
eigenvectors of the positive definite real symmetric matrix H (see e.g. [TQl Prop. 1.3.4]), 
shows that 

(Pll\GG T )} GeR02n ^ TlGG T /2 ^ = (P^\x 1 , . . . ,S2n))ME li2 „(^ 1 /2 c -/2) 

= P 1 ( 1 2) (l 2n )(2n), (3.6) 

where the final equality follows from (12.81) . Substituting this in (I3.5p . then substituting the 
result in the second formula of (11.41) shows that Q2n+i{ z ) = z 2n+1 — 2nz 2n ~ 1 as stated in 

m- 

We now turn our attention to the quaternion real Ginibre ensemble. Substituting (12. 1 D 
in the first formula of (11.51) gives 

2n 

QUZ) = E^(- 1 ) / ^( G ))G eQ S 2 „(e— 1 /2) - 

z=o 

Use of (12. 5p tells us that the average is non-zero when I is even only. Substituting this 
non-zero value gives 

n z 2(n-j) 

Qln{ Z ) = p(l/2)/ 1n x 2 ( GGr ^G€QS 2n (e-^GGt/2y (3.7) 

j=0 U ) 

Changing variables to H = GG\ which introduces the Jacobian factor det if (see [TQl 
Prop. 3.2.7]), then changing variables to the eigenvalues and eigenvectors of the positive 
definite self dual quaternion real matrix H shows that 

{P\J ' ) ( ( ^ ( ^ t ))GeQS 2 „(e- TrGGt /2) — {P\J ■ ■ ■ , x n))ME i>n (xe-*) 

= j p(l/ 2 ) (1 n ) -Q (2(n _ /+1))) (gg) 
1=1 



6 



where the final equality follows from (12. 8p . Substituting in (13. 7p . we find after simple 
manipulation that the second formula of (13.21) is reclaimed. 

With regards to the first formula in (13. 2p . ( 13 .4|) with R02 n (e~ TrGGT / 2 ) replaced by 
QS2n( e ~ Tl ' GGt//2 ) , we see from (12. 5p that only the average over Sii+i(G), I + 1 even, gives a 
non-zero contribution. Thus we have 

n 2(n-j) 

(det(zI Ar -G)TrG> G6QS2n(e ^ GGt/2) = -^^^^^^Ge^Me-^^y (3.9) 

Substituting (13. 7p and (13. 9p in the second formula of (11.51) shows cancellation occurs, leaving 
us with Q^ n+1 (z) = z 2n+1 as stated in (13721 . 

3.2 The induced Ginibre ensemble 

Let G r MxN , (M > N) denote the M x N rectangular analogue of the construction of a 
real Ginibre matrix. With O £ O(N), it is known [5] that the random matrix G] : = 
0((G\ IxN ) T G\ IxN ) l l 2 has a joint PDF proportional to 

(det(G^) T G^)( M -^/ 2 e- TrG ^ G ^ T / 2 . (3.10) 

A similar result holds true for quaternion real Gaussian matrices where, with 5* £ Sp(2iV), 
one has that [T^] the random matrix G} := 0((G r MxN ) T G T MxN ) 1 ^ 2 has a joint PDF pro- 
portional to 

(det(G?) T G?) M " 7V e- TrG ?( G ?) t / 2 . (3.11) 

Thus in both cases the joint PDF for the classical Ginibre ensembles is generalized by 
the additional factors (det GG^) a , for appropriate a. We can use (II. 4p and (II. 5p . now 
with G chosen from the corresponding induced ensemble, to calculate the skew orthogonal 
polynomials. 

Proposition 2. With a = (M — N)/2, we have for the induced real Ginibre ensemble 

Ql n (z) = z 2n , Ql n+1 (z) = z 2n+1 - 2(n + a)z 2n ~\ (3.12) 
and with a = M — n, we have for the induced quaternion real Ginibre ensemble 

n 21 

Ql + i(z) = z 2n+ \ Ql(z) = 2 n F(n + 1 + a) £ ^ (3.13) 

z=o ^ ' 

Proof. The derivation of the first result in ( 13 . 1 2[) goes through exactly as for the cor- 
responding result in (13. ip . which corresponds to the case a = 0. For the second result 
in (I3.12p . we see that (I3.5P again holds, provided we replace G £ i?0 2n (e~ TrGG I 2 ) by 
G £ R02n((det GG T ) a e~ TrGG I 2 ). Doing this implies that (13. 6p must be correspondingly 
modified, 

(P 1 ( 1 2) (GG r )) GeJ j 02ti(detGGT)0e _ TrGG T /2) = (P^\x 1 , . . . ,X 2 n))ME 1 , 2n (x- 1 / 2 +-e— / 2 ) 

= P 1 ( 1 2) (l 2 ")(2(n + a)), (3.14) 
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where the final equality follows from (12. 8p . Substituting in the second formula of ( II. 5p 
gives the second formula in (I3.12p . 

Turning our attention now to the induced quaternion real ensemble, we see that ( 13. 7p 
again holds, but with G G QS 2n (e" TrGGt / 2 ) by G G QS 2n ((det GG^) a e~ TrGG ^ 2 ), and that 

(P[J ) ( ( ^ ( ^ t ))G'6QS 2n ((detGGt) Q e- TrGGt /2) = {P\J 'fer • • .^DmE^Ii 1 ^- 1 ) 

= p(l/2) (1 n ) -Q (2(rl _ / + 1+a)); (315) 
1=1 

where the final equality follows from ( 12. 8p . Substituting in (13. Tj) gives the second formula 
in flffi]). 

The derivation of the evaluation of Q2 n +i( z ) proceeds in exactly the same way as 
done above for the evaluation of this polynomial in the classical quaternion real Gini- 
bre ensemble. Again there is cancellation, and we are left with the monomial evaluation 



4 Spherical and anti-spherical ensembles 

The strategy used to derive the evaluations of the skew orthogonal polynomials for the 
Ginibre and induced Ginibre ensembles can also be carried out for the real and quaternion 
real versions of the spherical and anti-spherical ensembles as specified by the joint PDFs 
(II. ip and (II. 2p . We consider first the case that G is real. 

Proposition 3. For the joint PDF U.l\) with G real we have 

Qln(z) = A Ql n+1 (z) = Z 2U+1 ~ ^ "\ + \ : /t?y (4.1) 

a 2 - (2n + «i + 1/2) 

while for the joint PDF U.S\) with G real we have 

%.(*) = *», M^-^ ji^i, , («) 

Proof. The mechanism for the evaluation Q\ n {z) = z 2n in both cases is precisely the 
same as our derivation of the first result in (13.11) . 

Consider then the evaluation of Qln+ii 2 )- The equations (13. 4p and (13. 5p again hold but 
with R0 2n (e" TrGGT/2 ) replaced by 

R0 2n ((detG t G) Ql /det(I + G t G') a2 ) and R0 2 „((det G ] G) h det(I - G ] Gf 2 ) 

for the spherical and anti-spherical ensembles respectively. We similarly modify (13. 6p so 
that in relation to (11. ip it reads 

(P 1 ( 1 ) (GG' T )) GgjR02n ( dctG tG) a i/det(I+Gl'G) O! 2) = (P[i{X\, . . . , X2n))ME 1/2 n{x- 1 ^ +a ^(l+x)- a 2) 

= Pl?(l 2n ) -r-— — — --, (4.3) 

11 V ; a 2 -(2n + «! + l/2)' V ; 
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where to obtain the final equality use has been made of ( 12.91) . and that in relation to ( II. 2p 
it reads 



(P X J \GG T )) GeR02n ( dctG >, G ) i a ldet ( I _ G <, G ^ 2 — [P[J \xi, . . . , X2n))ME 1 ,2n{x- 1 / 2+ ^(l-x)P2) 

11 1 j /3 2 + (2n + /3 1 + l/2)' l4 ' 4j 

where here the final equality follows from ( 12. 7ft . Substituting these formulas in the analogue 
of ( 13.51) . then substituting the result in the second formula of ( 11. 4ft gives the stated formulas 
forQ r 2n+1 (^). □ 
The real spherical ensemble, in the case a\ = 0, a 2 = N was first studied in [T2j . 
However there a fractional linear transformation was used to map the domain from the 
real line and upper half plane to the unit circle and unit disk. A direct comparison with the 
results (14.11) is therefore not possible; in fact the corresponding skew orthogonal polynomials 
were found to be monomials. However, in the recent Ph.D. thesis [I] a study of the induced 
real spherical ensemble has been undertaken. Then there is no advantage in introducing a 
linear fractional transformation. By making use of working which has a number of steps in 
common to that presented above, it is found that in the case of the parameters ol\ = L/2, 
a 2 = {m + L + N)/2 [3J Eqns. (4.2.97), (4.2.98)] 

Q^(x) = a* Ql j+1 (x) = x^ - Jl^j-l 3 ^' 1 - (45) 

This is precisely what is given by (14. ip with these choices of ai,a 2 , and N t— > 2n in the 
latter. 

The anti-spherical ensemble, with fti = 0, (3 2 = {L — N — l)/2 has been studied in 
[2T1 ITU1 [23] . and in [9] the skew orthogonal polynomials were evaluated as 

Q\ 3 {x) = x*\ Ql J+l (x) = x 2 ^ 1 - j^Tj'-' ■. (4.6) 

This is indeed the same result as is obtained by setting (3\ = 0, (3 2 = (L — 2n — l)/2 in 
( 14. 2ft . The corresponding induced ensemble has been studied in pE], corresponding to the 
parameters f3\ = L\/2 and (3 2 = (L 2 — N — l)/2. Indeed, substituting these values in ( 14. 2p . 
the latter with iV i— >■ 2n, we reclaim the explicit form of the skew orthogonal polynomials 
found therein [1] [(4.2.95), (4.2.96)]. 

We conclude with the calculation of the skew orthogonal polynomials for the spherical 
and anti-spherical ensembles in the case that G is quaternion real. 

Proposition 4. For the joint PDF U.l\) with G quaternion real we have 

r(2n + | + at - a 2 ) T(j + 1 + a<i) 

(4.7) 
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while for the joint PDF U.2\) with G quaternion real we have 

Proof. The mechanism for the result Q^n+iW = z 2n+1 in both cases is precisely the 
same as that already revealed for the corresponding result in (13.21) . 

For the evaluation of Q2 n { z )> the expression (13 .7p again holds, upon replacing the 
quaternion real Ginibre ensemble QS 2n (e~ TrGtG / 2 ) by the ensemble of 2n x In matrices 
with joint PDFs (11. ip and (II. 2p as appropriate. Furthermore 

„(det GtG) Q i / det(I+GtG) Q 2 ) 
= {P\J • • • i X n))ME iin {x 1+2a i{l+x)~ 2a 2) 

= ps /2 \n(-iy r ! Ql + t + m2n -] +1 * +ai - a j (4.9) 

13 V A ; r(2n + i + ai -a 2 )r(n-j + l + « 1 )' V 7 

where to obtain the final equality use has been made of (12. 9p . and 

(P 1 ( J 1/2) (G , G T )) G6 Q S2n ( dctG tG)/ 3 i det(I-GtG)^2) 
= (P^J (Xi, . . . ,a:n))ME 4 ,„(x 1 + 2 ' 3 l(l-a;) 2 ' 3 2) 



p (i/2) (r ) m + n + l)r(2n - j + \ + ft + ft) ^ 



1J ^ J T(2n + | + ft + ft)r(n — j + 1 + ft 
with the final equality now following from ( 127 
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